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ABSTRACT 

Using a high-resolution cosmological hydrodynamic simulation, we present a method to constrain 
extinction due to intergalactic grey dust based on the observed magnitudes of distant Type IA supernovae. 
We apply several simple prescriptions to relate the intergalactic dust density to the gas density in the 
simulation, thereby obtaining dust extinctions that may be directly compared to the observed distribution 
of supernova magnitudes. Our analysis is sensitive to the spatial distribution of grey dust, but is not 
dependent on its intrinsic properties such as its opacity or grain size. We present an application of our 
technique to the supernova data of Perlmutter et al., who find that their high redshift sample is ~ 0.2 
magnitudes fainter than the expectation for a non-accelerating, low-density universe. We find that for 
grey dust to be responsible, it must be distributed quite smoothly, e.g., tracing intergalactic gas. More 
realistic dust distributions, such as dust tracing the metal density, are inconsistent with observations at 
the 1.5 — 2cr level. Upcoming observations and improved modelling of the dust distribution should lead 
to stronger constraints on intergalactic grey dust extinction. 

Subject headings: Cosmology: observations, dust, extinction, large scale structure of Universe 



1. INTRODUCTION 

Recent observations o f Type IA Supern ovae (SNe) at 
redshifts up to z ~ 0.8 (Ricss ct al. 1998; Perlmutter et 
al. 1999, hereafter P99] ) have made possible classical cos- 
mological tests that require standard candles, such as the 
magnitude-redshift relation. The most dramatic result is 
that these SNe appear dimmer (by ~ 0.2 magnitudes) at 
high redshift than would be predicted in a non-accelerating 
universe, suggesting at face value that we live in an acceler- 
ating universe. However, other explanations are possible, 
including the one we consider here, namely that distant 
SNe are dimmer due to extinction by intergalactic dust. 

As distant standard candles, SNe are sensitive probes 
of extinction in the intergalactic medium (IGM). The dis- 
tribution of matter in the IGM may now be modelled ac 
curately in the context of moder n cosmology using hy- 



drodynamic simu l ations (see e.g., pen ct al. 1994 : Hern 
quist et al. 1996| ; pave et al. 1999|). The resulting IGM 



is not smooth, but rather traces large-scale structure. If 
such structure contains not only dark matter and gas but 
also dust, this would result in significant variations in SNe 
brightnesses due to intervening extinction. Observation- 
ally, the distribut ion o f Type la SNe magnitudes has a very 
small dispersion (P99). Thus by comparing simulations to 
the distribution of observed brightnesses, we can set limits 
on the amount of dust extinction and possibly constrain 
its spatial distribution with respect to intergalactic gas. 
In this Letter we present a techniq ue fo r doing this, and 
apply it to the SNe observations of [P99| . 

Intergalactic grey dust has been exa mined in a scries of 



papers by Aguirre (1999a, b; hereafter A99| ) and Aguirrc] 
fc Haiman (1999)| , who develop a scenario in which small 
grains are preferentially destroyed during ejection from 
galaxies, polluting the IGM with large dust grains that are 



effectively grey in the bandpasses of the SNe data. This 
greyness is necessary in order not to violate tight limits on 
reddening from SNe data, which imply th at ga lactic-type 
dust would provide negligible absorption (P99). Further- 
more, a significant fraction of the dust must reside in the 
IGM. If the grey dust causing extinction were present only 
in the ISM of the supernova host galaxy, this would intro- 
duce too large a dispersion in observed SNe magnitudes 
(Riess et al. 1998). In this study, we assume that grey dust 
blends smoothly from galaxies into the surrounding IGM. 
Our analysis is insensitive to the intrinsic properties of the 
dust, such as its opacity and grain size, since we use the 
simulations to directly translate the observed SNe mag- 
nitude distribution into a dust extinction in magnitudes. 
It is, however, sensitive to the way in which dust traces 
the distribution of gas in the IGM, and we will consider 
several simple but plausible variations of this relation. 

In § U we describe our simulations of the IGM, and of 
dust extinction. In § ^| we describe our analysis method 
and results, including constraints on grey dust afforded 
by current SN observations. In § |^ we discuss systematic 
uncertainties, and the implications of our results. 

2. SIMULATIONS OF INTERGALACTIC DUST 

We employ a hydrodynamic simulation of a A- 
dominated cold dark matter model, with VL m — 0.4, 
n A = 0.6, n b = 0.02/i- 2 , H = GSkms- 1 Mpc -1 , and 
as = 0.8. Our simulation volume is 50ft- -1 Mpc with 
lOh^ 1 kpc spatial resolution, having 144 3 dark matter and 
144 3 gas particles, and was evolved from z = 49 — > u sing 
Parallel TreeSPH flDave, Dubinski fc Hernquist 1997[ ). 

In order to obtain dust column densities along lines of 
sight, we consider three different ways that dust may trace 
intergalactic hydrogen gas: 
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1. Pdust oc p gas : Dust traces gas linearly. 

2. pdust oc Pmctai: Dust traces metals linearly. 
3- Pdust oc Pg as : Dust traces gas quadratically. 

As our simulation makes no direct prediction for the metal- 



the constraints that can be put on the dust extinction by 
using the observed distribution of SNe magnitudes. 



licity of gas, w e adopt a heuristic prescription (c.f. Cen & 
Ostriker 1999| ), in the second case above. We assume that 



the metallicity is 10 2 solar if the gas overdensity is less 
that 10, solar if the overdensity is greater than 1000, and 



g-lincar in between. 


( a ) Pdust^gas 




- 




z = 0.05 


z = 0.5 



( b ) Pdust "^gas 
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z = 0.05 





Figure 1: Dust extinction maps in 2.2° x 2.2° patches of sky out 
to z = 0.05 (left) and z — 0.5 (right), for (a) pdust oc p gaa and 
(b) pdust oc pg aa The median extinction magnitude to z = 0.5 
is set equal to 0.4 for both (a) and (b). The maps correspond 
to what would be seen against a white background. 

To extract dust extinction values from the simulations, 
we assume that the gas associated with each part icle is 



spread over its SP H smoothing volume (see e.g., Hern 
quist & Katz 1989). We perform a numerical integration 
of gas column density along 5000 rays cast through these 
volumes, at the same time applying one of the three trans- 
formations given above to relate gas to dust densities. To 
reach the required path lengths out to z ~ 0.5, we fol- 
low rays through 26 simulation volumes, each ray entering 
through a random point on a random face. This yields the 
column density of dust to z = 0.5 along each line of sight. 

In Figure 1 we show extinction maps of 2.2° x 2.2° 
patches of sky, for (a) p dus t oc p gas and (b) p dus t oc p 2 as . 
The median extinction to z = 0.5 was set to be equal 
(to 0.4 mag) for cases (a) and (b). In Figure 2, we show 
how the mean dust extinction and its dispersion varies 
with redshift. Here, the dispersion in all three panels was 
set to be the same small value, equal to the difference in 
quadrature of the dispersion in SNe magnitudes at high 
(o"z=o.5 = 0-15 7) and low (c z= o.o5 — 0.154) redshifts ob- 
served by P99, namely 0.03 mag. Figs. 1 and 2 show that 
the mean extinction is much greater when the dust is more 
smoothly distributed. We will now quantitatively explore 

1 We use the median in order to be less sensitive to long tails of the 
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Figure 2: Mean and dispersion of the dust extinction as a func- 
tion of redshift. Filled squares show mean, open circles with 
line show the dispersion. Values have been scaled to a disper- 
sion of 0.03 magnitudes at z = 0.5. 



3. COMPARISON WITH OBSERVATIONS 

We make use of two characteristics of the observed SN 
data in our comparison, the change with redshift of the dis- 
persion in SNe magnitudes, and the shape of th e histogram 
of SNe magnitudes. As mentioned above, P99 found little 

0.05 



As mentioned above 
difference in the dispersions of two samples with z 
and z ~ 0.5. As they stated, this leaves little room for 
dispersion due to dust, as this dispersion is expected to 
increase for longer path lengths (see Fig. 2). In order to 
quantify this, and the effect of the distribution shape, we 
gen erate simulated SNe magnitudes, and compare them to 
the |P99] data using a maximum likelihood approach. 

T he observational datasets we use are both taken from 
P9| (their Tables 1 and 2), being the high-z SNe of the 
Supernova Cosmology Pr oject, and the low- z sample of 
Calan-Tololo SNe survey flHamuy ct al. 1996| ). We use 40 
(non-reddened SNe) of the former SNe, between z = 0.172 
and z = 0.83 (z ~ 0.5) , and 16 of the latter SNe which lie 
between z = 0.02 and z = 0.101 (z ~ 0.05) . 

We generate simulated datasets for each dust model de- 
scribed in § 2. For each dust model, we vary two parame- 
ters; first, Mc, a cosmological magnitude shift applied to 
all simulated SNe at a given z, normalized so that Mc = 
at z = .5 corresponds to the best fitting cosmology found 
by [P99] (with Q, m = 0.28 and f7 A = 0.72); and second, A v , 
the median V-band magnitude of dust extinction out to 
z — 0.5 1 . Mc ~ —0.2 then corresponds to an open model 
with TO ~ 0.3, and Mc ~ —0.4 to an TO = 1 model. We 
generate the simulated datasets as follows: 

(1) We renormalize the 5000 simulated lines-of-sight so 
that the median dust extinction to z — 0.5 equals Ay- 

(2) We add the cosmological shift, Mc- 

(3) We broaden the magnitude distribution, which involves 
convolving it separately with a Gaussian of width given by 
the observational error of each SN, and then adding these 

distribution (c.f. Figure 3) when renormalizing extinction values. 
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distributions together. When doing this, w e inc lude an 
"intrinsic" SN dispersion of cr int =0.17 mag ( P99| ). Vary- 
ing this (by ±0.1 mag) makes little difference to the results 
(see §4). 

(4) We truncate the distribution at a magnitude difference 
Am = 1, to roughly account for the fact that SNe along 
lines of sight passing through high extinction regions would 
not make it into these magnitude-limited samples. 




Figure 3: The PDF of magnitude differences, Am, for 5000 
simulated SNe. The shaded region is included in the calcu- 
lation of the dispersion; non-shaded region indicates heavily 
obscured SNe that are not included in the PD F no rmalization. 
The histogram shows 40 of the high-z SNe of P99. 



We derive the probability distribution function (PDF) 
of SNe magnitude differences Am predicted by the simu- 
lation, so that the predicted number of SNe between Am 
and Am + cIAm is NP(Am)cIAm, where N is the number 
of observed SNe. Figure 3 shows histograms of the PDF 
of Am, for Ay = 0.4 ma g, and Ale = —0.4. We also plot 
the observational data of P99 . In effect, for this plot, we 
have brightened the simulated SNe to mimic an fi TO = 1 
model and then dimmed them with dust. We can see that 
in panel (a), where the dust is fairly smoothly distributed, 
the simulated PDF is not too different from the observa- 
tions. In the other panels, which have clumpier dust, there 
is a skewness not seen in the observed data, due to a long 
tail of high extinction lines-of-sight. 

For each set of simulated lines of sight, we form the 
relative likelihood of drawing all the observed SNe mag- 
nitudes, C — P(Am,i), where Ajj,i is the magnitude 
difference of SN i. We define the quantity S = — 2 ln£ and 
assume that S follows a % 2 distribution in order to derive 
confidence limits on the parameters Mc and A v . In the 
present analysis, we use results at two different redshifts 
(z = 0.05 and z — 0.5), and combine the two by adding 
the values of S. It is this step, combining the likelihoods 
at two different redshifts, which constrains the amount of 
additional dispersion (or the change in the shape of the 
magnitude distribution) due to grey dust between low and 
high redshifts. Contours of AS = 2.3,6.2, and 11.8 (the 
difference in S from its minimum value), representing 1, 2 
and 3cr intervals of joint confidence, are plotted in Figure 4. 
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Figure 4: 1, 2 and 3<r contours for the cosmological magni- 
tude shift, Ma, with respect to a model with Qa = 0.72 and 
Qm = 0.28, and the median extinction due to grey dust (Av)- 

Figure 4 shows that the smoothest distribution of dust 
easily accomodates enough extinction to reconcile a non- 
accelerating (Mq = —0.2) or flat (Mc = —0.4) universe, 
as there is a strong degeneracy between the cosmologi- 
cal magnitude shift and dust extinction. The other cases, 
where dust traces metals and /?g as , are coincidentally quite 
similar. Such a distribution of dust would be mildly disfa- 
vored in a non-accelerating open universe and ruled out at 
~ 99% confidence in an Einstein-de-Sitter universe. The 
constraints arise mostly due to the shape of the distribu- 
tion; the fact that the observed dispersions are similar at 
z = 0.05 and z = 0.5 is relatively unimportant, only mak- 
ing a noticeable difference in panel (a) where the shape is 
similar to the observed distribution. 

4. DISCUSSION 

We have found that current SNe datasets appear to have 
some power to constrain grey dust models. We find it 
somewhat unlikely (at the 1.5 — 2 level) that sufficient grey 
dust could be distributed in the relatively clumpy fashion 



expected fo r the metal distribution (c.f., Sncdin fc Os- 
trikcr 1997| , [Dave et al. 1998j , |Cen fc Ostriker 1999| ) to 
reconcile the |P99| data with a non- accelerating universe. 
This suggests that any substantial grey dust component 
must be largely segregated from the galax ies w here it was 
formed. In the specific grey dust model of A9£ , sputtering 
does destroy dust more effectively in denser regions, but 
once in the IGM the large grains have long lifetimes, so 
that dust is still likely to trace the metals. The question 
of how the dust and metals are distributed can be an- 
swered self-consistently by modelling the relevant physical 
processes (dust and metal ejection from galaxies) directly 
in the simulations (e.g., Aguirre et al., in preparation). 

There are many alternative explanations for SNe ap- 
pearing dimmer in the past. For exa mple, there may b e 
metallicity effects in the host galaxy ( Hdflich et al. 1999), 
intrinsic evolution (Riess et al. 1999), observational se- 
lection differences between CCDs used for distant SNe 
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and p hotographic plate s used for nearby samples (Howell 
Wang fc Wheeler 1999 ), or time evolution of the gravita- 



tional constant ([Amcndola, Corasaniti fc Occhionero 1999 



Garcia-Berro et al. 1999 ). Such effects would alter the in- 
terpretation of our paramete r Mc- Gravitational lensing 
magnification (Metcalf 1999) would also change the mag- 
nitude dispersion at higher rcdshifts. 

The study we have presented is reasonably general, in 
that our analysis is not dependent on the (unknown) mi- 
croscopic properties of this hypothetical intergalactic dust. 
Still, there are some possible systematic unc ertainties, 
which we now consider. We assume (as do P99| ) that the 
observed dispersion in SNe magnitudes in excess of the 
estimated observational errors is an intrinsic property of 
SNe and does not vary with redshift. One could envision 
scenarios in which the intrinsic dispersion is lower at high 
redshift, thus allowing more dispersion from dust. How- 
ever, as mentioned earlier, most of the statistical power 
of our analysis comes from the shape of the distributions, 
which is not significantly affected by a change in the in- 
trinsic dispersion. Also, we find that when we decrease the 
intrinsic SN dispersion (<Xi„t), models with more dust fit 
slightly better at high-z, but the low-z fit becomes worse, 
a trade-off which means that the overall results hardly 
change. The distribution shape depends on our assump- 
tion that the intrinsic dispersion has a Gaussian distri- 
bution in magnitudes, whereas a distribution skewed to 
fainter magnitudes could weaken our constraints. With a 
larger sample of low-z SNe we might be in a position to 
test this, by using the distribution of low-z SN magnitudes 
to make a simulated high-z sample with the correct intrin- 
sic distribution shape. Also, we have made a simplifying 
approximation in our simulated high-z samples, by using 
dust extinction magnitudes that result from integrating 
the dust contribution from z — to exactly z = 0.5. The 
real SNe are at a range of redshifts; our constraints are 
effectively conservative because of this. We have also as- 
sumed that the total mass of dust does not change with 
redshift from z — 0.5 to z — 0; dust increasing with time 
would strengthen our constraints, while a decrease seems 
implausible. Finally, the truncation of our simulated mag- 
nitude distribution at +1 mag is a rather approximate 
procedure. We have tried changing the cutoff, and find 
that with none (an unrealistic case) the constraints be- 
come stronger, as we might expect. With a lower cutoff, 
the right side of the observed distribution is not repro- 
duced. One additional degree of freedom would involve 



changing the functional form of the cutoff. In the future, 
as observations improve, we plan to simulate such obser- 
vational selection effects more carefully. 

On the simulation side, our relatively low mass resolu- 
tion (m gas = 8.5 x 10 8 M Q ) means that we miss fluctuations 
in the extinction that occur on smaller mass scales. If we 
had higher resolution, this should have the effect of making 
our dust constraints tighter. We find that the dispersion in 
projected magnitudes largely depends on small scale fluc- 
tuations. We tested this by splitting the simulation into 
small sub- volumes and then shuffling them to remo ve large 
scale correlations (see e.g., Bhavsar & Ling 1988) before 
projection, and we obtained similar results. Changing the 
assumed cosmological model could have some effects, al- 
though it is difficult to conceive of models which have less 
pow er on small scales wh ile fitting other constraints (see 



e.g., White fc Croft 200C) 



The models of |A99| have grains that do produce some 
reddening, typically in the infrared. A complementary ap- 
proach to ours is therefore to use information fro m di ffer- 
ent color bands. Such an approach was used by P99, for 
example, to constrain n ormal galactic-typ e dust. Recent 
near-IR observations by Riess et al. (2000) have made red- 
dening constraints tight even for non-standard large dust 
grains. Our approach does not make any use of color in- 
formation, and so constrains the most extreme scenario in 
which the dust is totally "grey" . 

In summary, we have used cosmological hydrodynamic 
simulations to explore how intergalactic grey dust could 
affect observations of high redshift supernovae, and how 
supernova data can constrain grey dust extinction. We 
conclude that only a fairly smooth distribution of dust 
could readily mimic the effect of an accelerating Universe. 
Such a distribution would be strange, as the dust would be 
strongly biased away from metal producing regions. More 
realistic dust distributions are mildly disfavored, but up- 
coming samples of SNe (at current redshifts) should enable 
us to put tighter limits as they more precisely determining 
the shape of the brightness distribution. At higher red- 
shifts (z ;> 1), grey dust predicts that SNe should show 
increased dimming, while decreased dimming would be 
strong evidence for a cosmological constant. 



We thank Anthony Aguirre, Bruce Draine, and Bob Kir- 
shner for useful discussions, and David Weinberg for help- 
ful comments on the manuscript. 
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